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bjo! Abstract 
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We study the nuclear stopping in high energy nuclear collisions using the 
constituent quark model. It is assumed that wounded nucleons with different 
number of interacted quarks hadronize in different ways. The probabilities 
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00 ■ of having such wounded nucleons are evaluated for proton-proton, proton- 

O 

nucleus and nucleus-nucleus collisions. After examining our model in proton 

ON i 

proton and proton- nucleus collisions and fixing the hadronization functions, 
q ■ it is extended to nucleus- nucleus collisions. It is used to calculate the rapidity 

distribution and the rapidity shift of final state protons in nucleus-nucleus 
collisions. The computed results are in good agreement with the experimental 

X 

^ . data on 32 S + 32 S at E lab = 200 AGeV and 208 Pb + 208 Pb at E lab = 160 

AGeV. Theoretical predictions are also given for proton rapidity distribution 

in 197 Au+ 197 Au at yfs = 200 AGeV (BNL-RHIC). We predict that the nearly 

baryon free region will appear in the midrapidity region and the rapidity shift 

is (Ay) = 2.22. 
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I. INTRODUCTION 



Whether the incident nucleons are stopped in or passed through the target nucleus is a 
fundamental and important concern in high energy heavy ion collisions. These two different 
situations emerge according to collision energy and atomic mass of nucleus. They suggest 
different relevant dynamics, the dynamics of shock formation and Landau hydrodynamics ]I[ 
in case of stopping regime or the Bjorken longitudinal expansion and inside-outside cascade 
dynamics [Ej in case of baryon free regime. The distinction between the two situations 
is related to how the incident nucleons slow down by multiple collision with nucleons of 
the other nucleus and how the collision energy are deposited for particle production. The 
relevant measure, the nuclear stopping power provides an estimation of the energy density 
achieved in collisions. Experimental data indicate so far a baryon rich regime at midrapidity 
for heavy nucleus collision such as gold on gold at AGS energies and sulphur on sulphur 
0] or lead on lead || at SPS, where the achievement of initial condition for QGP formation 
is still controversial. The clear baryon free regime may be realized at RHIC energy region. 

For the description of high energy nucleus collision dynamics, there are many theoretical 
models which differ in their assumptions as to how the particles share the incident energies, 
where the sources of particle production are and how the produced particles hadronize, while 
having almost the same picture for a multiple collision process of constituents such as the 
Glauber model ||. A detailed comparison of these models can be found in Ref. [|7| by Werner 
and Ref. |§ by Wong. In particular, for the proton rapidity density distribution which is 
related to how high baryon densities may be attained in reactions, these model predictions 
show notable difference due to their different assumptions. For example, as discussed by 
Gyulassy 0, the RQMD [|10J predicts a much higher degree of baryon stopping than VENUS 
|7j and HIJING |TT| which predict a concavity at midrapidity for a reaction such as lead 
on lead collision at SPS energies. In some approaches a nucleon-nucleon model is directly 
extended to nucleus-nucleus collisions without examining proton-nucleus collisions. Our 
opinion is that model for nuclear collisions should first be examined for NA collisions and 
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then extended to AB collisions. 

The notion of constituent quarks as the units of collision has been shown to be very 
useful to describe not only hadron-hadron collisions but also hadron-nucleus(/iv4) collisions 
at high energies. In particular, the successful results has been obtained in application to the 
projectile fragmentation region of HA collisions, since the first application to HA collisions 
by Anisovich et at fT2|] and, later on, further developed by several authors ||T3|-[T7| . 



In this paper, We first formulate the constituent quark model(CQM) relevant to nucleon- 
nucleon, nucleon-nucleus and nucleus-nucleus collision dynamics at high energy by consid- 
ering quarks, the constituents of hadrons, as the unit of collision instead of hadrons. When 
the projectile nucleon is incident on the target in nucleon-nucleon collision, the collided 
nucleon has the probability of becoming three types of wounded nucleon with one, two or 
all the three interacted quarks. In case of nucleon-nucleus collision, the incident nucleon 
after collision has these probabilities by multiple collision with individual nucleons within 
the target nucleus. In nucleus-nucleus collision, each collided nucleon from the projectile 
nucleus becomes one of the three types of wounded nucleons. This distinction of the three 
types of wounded nucleons is a crucial difference from multiple-collision models constructed 
at the nucleon level. The three probabilities of quark absorption can be calculated from the 
given nuclear density and the total inelastic quark-quark cross section of* eZ . 

The three probabilities are used to estimate the degree of nuclear stopping or the rapidity 
distribution in high energy nucleus-nucleus collisions. The rapidity density distribution or 
the momentum distribution in nuclear collision is expressed by three factors, the flavor factor, 
the quark interaction probability of the incident nucleon and the fragmentation function. 
The flavor factor is introduced to make a distinction whether a final state baryon is an 
observed proton, a neutron or a nucleon decayed from a hyperon. The determination of 
fragmentation function fi(x) comes from fitting the data on proton-proton ]TJ| and proton- 



nucleus collision [O. 



This paper is organized as follows. In section [TT], we show the assumption of CQM and 
the quark interaction probabilities having three types of wounded nucleons with one, two 



or all the three interacted quarks in nucleon-nucleon collision at high energy. In section |T|, 
the quark interaction probabilities in nucleon-nucleus collisions are given as a function of 
the mass number A in nucleon-nucleus collision. In section [IV], the average numbers of three 
type wounded nucleons are calculated for AA interaction. In section [V], we compute the 
fractional momentum distribution of protons in pA collision. In section [VT|, we calculate the 
rapidity distribution of leading protons in AA interaction. Finally, section [VI I| is devoted to 
conclusions and discussions. 
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II. QUARK INTERACTIONS IN NUCLEON-NUCLEON COLLISION 



We start giving the outline of the assumptions of the CQM pertinent to the description 
of nuclear stopping in high energy nuclear collisions. The CQM is based on three funda- 
mental assumptions which are related to the structure of hadrons, the interactions between 
constituents of the projectile and the target hadrons and the hadronization of quarks in 
participant nucleons JE2]-|llJ. We first assume that a hadron(meson or baryon) consists of 
two or three spatially separated constituent quarks. Secondly, in a hadron-hadron, hadron- 
nucleus or nucleus-nucleus collision, some quarks from the projectile are assumed to interact 
independently with some quarks from the target thus losing a considerable fraction of their 
initial momenta, while the quarks which escape from colliding in both the projectile and 
the target pass through retaining their initial momenta. The third assumption claims that 
those quarks hadronize eventually via fragmentation and recombination mechanism. 

In order to calculate the total inelastic cross section for nucleon-nucleon collisions in 
terms of CQM, we need the probability of having a quark-quark inelastic collision when one 
quark in the projectile is at an impact parameter b relative to another quark in the target 
which is given by 



in the point particle approximation. The integration over impact parameter gives the total 
inelastic cross section of quark-quark collisions; 



We consider the collision of a beam nucleon N with a target nucleon N'. To the proba- 
bility of a quark-quark collision Eq. ([I]), by multiplying the probability elements for finding 
a quark p N {b N )db N dz^ and pN>{b N i)db N idz N i in the volume element db N dzjy and db N /dz N >, 
and by integrating over the collision axis z^ and zn>, we obtain the probability that a par- 
ticular quark in N interact with a particular quark in N' when N and N' are at an impact 
parameter b relative to each other, 



(1) 




(2) 
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W(b) = J db N db N ,p N (b N )p N ,(b N ,)h(b - b N + b N ,), (3) 
where Pat(6at) is the ^-integrated normalized quark distribution in the nucleon; 

pN(b N ) = J p N (b N ,z N )dz N . (4) 

It is normalized as 

J p N (b N )db N = 1. (5) 

Using W(b) of (|]), we can evaluate the total inelastic cross section and various proba- 
bilities. When a projectile nucleon is at an impact parameter b relative to a target nucleon, 
the probability of the occurrence of an inelastic event is 

g(b) = l-{l-W(b)Y, (6) 

where the second term means the probability that all quarks in the projectile pass through 
the target nucleon without any inelastic collision. Therefore, we obtain the total inelastic 
cross section af n ^ for iVW collision 

*£2 = J dbg{b). (7) 

The probability g(b) of Eq. ([]) can also be expressed as a sum of g^(b) which is the 
probability that i quarks in the projectile nucleon iV collide with any quarks of the target 
nucleon N'; 

g(b) = j2g®(b), (8) 

8=1 

where 



gW(b) = n [1 - {1 - W(b)Y\" [1 - W(b)f 3 ^ . (9) 

Here (^j is the combinatorial factor and the term 1 — {1 — W(b) } 3 implies the probability 
that a particular quark out of the projectile nucleon interact at least once with the quarks 
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in the target nucleoli. Integrating over the impact parameter b and dividing by the total 
inelastic cross section, we obtain three integrated probabilities that i quarks out of three are 
absorbed in nucleon-nucleon collision; 

«... = J— ( Aha® 



P { n>n> = ^W dbgW(b), (10) 

a inel J 

which is illustrated in Fig. [I]. 

Given the quark density distribution pjv(6) and the inelastic quark-quark cross section 
cr^ei, we can calculate the total inelastic proton-proton cross section a^ el and the probability 
of quark absorption Pj£j in proton-proton collision. The quark distribution is assumed to be 
Gaussian for simplicity, 

M&) = 2^exp(-A_), (11) 



where the parameter (3 is related to the root-mean-square radius of the proton r£ 



rms ' 



P 2 = l^rmsf- (12) 

Electron scattering data [^] gives r v rms = 0.862 fm. For the inelastic quark-quark cross 
section as the input, we take af^ el = 4.32 mb in order to reproduce the total inelastic cross 
section of pp collision <jf^ el = 30 mb which is observed for center-of-mass energy 3 GeV ^$ 



< 100 GeV and a™ el = 6.64 mb to give o% = 41 mb for y/s = 200 GeV f2|. The 



numerical values of the probability of quark interaction in pp collision for cr^ei = 4.32mb are 
shown in Table |. It should be noticed that Pj® is considerably large implying violation of 



the additive quark approximation [22 . 
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III. QUARK INTERACTIONS IN NUCLEON-NUCLEUS COLLISION 



When an incident nucleon collide with the target nucleus, the projectile nucleon can 
interact with many nucleons in the nucleus. We use the probabilities of quark absorption in 
NN collision to obtain those in NA interactions. When a nucleon is incident at an impact 
parameter b relative to the nucleus A, the probability for the nucleon to collide with a 
particular nucleon in the target nucleus is given by 

V A (b) = J db A p A (b A )g(b - b A ), (13) 

where p A {b A ) is the ^-integrated nucleon density distribution of nucleus A, db A p A {b A ) is the 
^-integrated probability element of finding the nucleon and g{b — b A ) is the probability of 
inelastic NN interactions at impact parameter b — b A given by Eq. (§). 
The total inelastic cross section for NA collisions is given by 



NA 
inel 



db 



l-{l-V A (b)} A \. (14) 
As g(b — b A ) in Eq. ([13]) is given by (||), V A {b) is decomposed into a sum of three terms; 

W)=Evf(b) (15) 



i=l 



where V\ (b) is the probability that i quarks in the projectile nucleon interact with a nucleon 
in the target nucleus, 

vf(6) = J db APA (b A )g®(b - b A ). (16) 

Let us calculate the probabilities Pj^ A for projectile nucleon having j interacted quarks in 
NA collisions. We first expand of Eq. (14) into a sum of the contributions from n 
N — N collisions; 

°22 = E / db(fj [v A {b)f [i - v A (b)] A - n . (ir) 

By substituting Eq. ([TBI) f° r (PD and expanding the latter, we obtain 



S 



p, 



0) 



E 



'A) 



°inel n=l \ n > 



dbUtf A (n;b){l-V A (b)} 



A-n 



NA — ^NA 

where U^ A (n; b) is the probability of having j interacted quarks in n N 
is given by 

U$ A (n;b) = U§(n;b); 

U% ] A (n-,b) = U$(n;b) + U$(n;b) 



(18) 

N collisions and 

(19a) 
(19b) 



n-l 

+ E 

k=l 



n 

. k , 



b)U$(n -k;b) + ±U$(k; b)U$(n - k; b) 



4 3 i(n; b) = [V A (b)] n - U$> A (m b) - U% A (n; b) 



(i) 



'(2) 

r(2) 



where U$ (n; 6) is the probability of having j interacted quarks in n N — 
having i interacted quarks in each collision and is given by 



(19c) 

^collisions while 



C/W(n;6) = 3(-y«(6) 



2)1^(6) 



C/g}(n;6) = 3(M 2) (&) 



(20a) 
(20b) 
(20c) 



r ( ( §(n;b) = 3(2 n - ' 
(2) (n;6J = 3 (J 1 

In the Appendix [S] we present the details of our calculation of P^\. 

Now we need to fix the nucleon density distribution p A {b) in order to carry out the 
numerical calculation. We use the Woods-Saxon parameterization, 

Po 



pA{r) 



exp 



r-R\ 



for heavy nuclei (A > 16) while the Gaussian distribution 



Pa{t) 



1 



(2vr/3 2 ) 3 /2 



for 9 Be, where both distributions are normalized as 



d 3 rp A (r) 



(21) 



(22) 



(23) 



The parameters for each distribution can be obtained from the data on elastic electron and 
hadron scattering on nuclei. For the Woods-Saxon distribution, we use the parameter R 



and a in Ref. 23 : 



^=1.12^3-0.86^-3, a = 0.54fm (24) 
for (A > 16). For the Gaussian distribution, we take [24]: 



P 2 = \{rf ms )\ = 2.3fm (25) 

for 9 Be. The calculated values of the total inelastic cross section af^ei an d the probabilities 
PpA, Pj$ and are shown in Table | and the behavior of P^jj as a function of A is shown 
in Fig. |||. The calculated cross sections are in good agreement with the experimental values 



1^ , 25] . The process of one quark interaction with the probability pty is dominant for light 



pA 

nuclei. However, even for Be, the probability of having two interacting quarks is not so 
small. All the three probabilities have comparable magnitudes for A > 60. 
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IV. QUARK INTERACTIONS IN NUCLEUS-NUCLEUS COLLISIONS 



In nucleus-nucleus collision, many nucleons of the projectile nucleus(A) may collide with 
many nucleons in the target nucleus(-B) by multiple collisions. After the full process of 
multiple collisions, the incident nucleons are divided into four types of nucleons having zero, 
one, two or all the three interacted quarks. When the beam nucleus A and the target 
nucleus B are situated at an impact parameter b relative to each other, the probability for 
the occurrence of a collision between one nucleon in A and one nucleon in B is given by 

V AB (b) = J db A db BPA (b A )p B (b B )g(b + b A - b B ), (26) 

where g(b) is given by Eq. ([]). 

The total inelastic cross section for AB collision is given by 



°t3i = I db 



inel 



1 - {1 - V AB (b)} AB . (27) 



In the same way as in Eqs. (|15|) and fllED , the probability V AB {b) can be expressed as a sum 
of three probabilities of incident nucleons having i interacted quarks 

3 



v AB {b) = Y,nm (28) 

i=l 

where 

Vffi(b) = / db A db BPA {b A )p B {b B )gV(b + b A - b B ), (29) 

with g^\b + b A — b B ) being given by Eq. (|9|). The expansion of total inelastic cross section 
Eq. ( P?D gives the probabilities of having m wounded nucleons on the beam side 

PabM = 4f / db ( A ) f 1 - i 1 - VAB(b)} A ] m [1 - V AB (b)f^ m) , (30) 

a inel J \ m J 

and the average number is 

A 

(m) = rnP AB (m). (31) 



m=l 



The number m is a sum of mi, the number of nucleons having i interacted quarks; 
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m 



i=l 



(32) 



We shall evaluate the probability Pab{™>\ TOi, TO2, TO3) of having such configuration of 
wounded nucleons. It is convenient to introduce the probability R^ B {b) that one of in- 
cident nucleons has j interacted quarks at a fixed impact parameter b. It is obtained in a 
similar way as in NA collision by expanding the factor 1 — {1 — V AB {b)} A . The result is 



with 



Rf B (b) = £ [ B W(n; b) [1 - Vis(6)}] 

n=l \ n / 



Y.R { l ) B (b) = i-{i-v AB (b)} B . 

3=1 



B-n 



(33) 



(34) 



Here, U^ B (n; b) is obtained by substituting V AB (b) in AB interaction instead of (b) in 
NA interaction in Eqs. (|T9"|) and (p0|). Using the polynomial expansion, we have 



P AB (m;m 1 ,m 2 ,m 3 ,) 



AB 

u inel 



db 



R%iP) - R%{b) R%{b) {1 - V AB (b)} 



'A" 
m / 

mi 



m 

>(2) 



m — mi 

y m 2 , 



in-2 



(3) 



B 



A-m 



(35) 



where m 3 = m — mi — m 2 . The average number of mj is given by, 

A m m—mi 

( m j) = J2 J2 m j P AB (m;mi,m 2 ,m 3 ). (36) 

m=l mi=0 m2=0 

Using the nucleon density distribution of Eq. (|2"ID with Eq. (|24j) , the nucleus- nucleus 
inelastic cross section and the average number of wounded nucleons are calculated for 32 S 
+ 32 5 ( Elab = 200 AGeV), 208 Pb + 208 Pb {E lab = 160 AGeV) and 197 Au + 197 Au (^ = 200 
AGeV) reactions as there are recent experimental data from CERN SPS for the first two 
reactions. The third reaction is chosen to give a prediction for RHIC experiment. Along 
with a{^ A h the calculated values of (to), (toi), (to 2 ) and (771.3) for no-bias events are listed 
in Table [TJ. For Au + Au and Pb + Pb reactions, we have shown the values estimated 
from the empirical formula (J^nei = nr^A* + B~$ — b) 2 with parameters r = 1.48 fm and 
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b = 1.32 fm [26]. The calculated cross section is larger than the experimental values for 



S + S, while it is in good agreement with data for Au + Au and Pb + Pb reactions. The 
discrepancy of the cross section in S + S does not affect on Pab^] mi, ^2, ^3) and (rrij). 

The probability distributions Pab(^) of Eq. ( |30| ) and the corresponding average values 
(m) for S + S and P6 + Pb reactions for different triggers are shown, respectively, in Fig. 
H] and Table || where the values of (uij) are shown also. Here, Veto and Et triggers in S + 
S reactions mean central triggers which correspond to 2% and 11% of total inelastic cross 
section, respectively [|J. Corresponding to these triggers, we introduce the impact parameter 
cutoff b max in our theoretical calculations so that the partial cross section 

/ dbaf n B el (b) (37) 

J \b\<b max 

becomes 2% or 11% of the total cross section. The cutoff is found to be b max = 1.18 fm for 
the Veto trigger and b max = 2.78 fm for the Et trigger. For both Pb + Pb and Au + Au 
reactions, we take b max = 6.05 fm and b max = 5.86 fm, respectively, corresponding to 15% 
of of^ v On the other hand, the peripheral trigger in S + S reactions correspond to 65% of 
a tneii which lead to the small \b\ cutoff b min = 4.9 fm. 

In Fig. |3](a), the distribution P^_b(^) for no-bias event of S + S reactions shows the 
well-known horse-back shape. It is obvious that the distribution at small m is dominated 
by peripheral collision while the one at large m is by the central ones. The clear separation 
of two components in the m-distribution is the consequence from the strong correlation 
between m and b. As shown in Fig. 0(b), in Pb + Pb collisions, the distribution for no-bias 
events show prominent peaks in both the smallest and the largest m-region. The largest 
m-region is of course dominated by the central collisions. The most remarkable features of 
the corresponding average value (m) shown in Table [IT] is that more than 94% of the incident 
nucleons are wounded in the central collisions of both Pb + Pb and Au + Au. 

Shown in Fig. ^ are m-dependence of fractions (rrii(m)) /m in no-bias event of S + S 
and Pb + Pb reactions. Here (mi(m))/m is the average value of for a fixed m. At small 
m, the wounded nucleon having one interacted quark mi dominate in both reactions. At 
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larger m, the fractions (mi(m))/m and (m2(m)}/m are comparable to each other in S + S 
reaction. On the other hand, for m > 197, the fraction {m^{m))/m becomes largest in Pb 
+ Pb reactions. For this region, (m.3) is largest in central Pb + Pb collision at y/s = 17.4 
AGeV as seen in Table |IT[ The trend is stronger in Au + Au collisions at RHIC energy 
because of larger cr^ei- 
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V. FRACTIONAL MOMENTUM DISTRIBUTIONS OF PROTONS IN 
PROTON-NUCLEUS COLLISIONS 



In this section, we study the inclusive spectra of protons in the projectile fragmentation 
region of proton-nucleus collisions as it gives a basic information on the nuclear stopping 
in high energy nuclear collisions. According to CQM, both participant quarks and specta- 
tor quarks hadronize via fragmentation and recombination. However, instead of giving a 
detailed description of such hadronization dynamics, we here introduce a phenomenological 
fragmentation function fi(x) (i = 1, 2, 3) which is fractional momentum distribution of pro- 
tons coming from a wounded nucleon having % interacted quarks. Assuming the independent 
fragmentation of each type of wounded nucleons, we express the proton spectra in pp, pA 
and AB interactions as 

3 



dN 
dx 

dN 
dx 

dN 
dx 



= EV^ (38a) 

pp— >pX i=l 

= EV$/<(*), (38b) 

pA— *pX i=l 

= j2\^(m t )f l (x), (38c) 

AB^pX i=l 



for < x < 1, where dN/dx is the normalized single particle inclusive cross section 
cruder/ dx and x is the Feynman scaling variable defined in c.m.s. Moreover, Aj and A^ 
are the flavor factors which can be interpreted as the probabilities of finding a proton in the 
hadronization product from wounded nucleon having % interacted quarks provided the effect 
of baryon-antibaryon pair production is negligible. The probabilities and P^\ are given 
by Eqs. ([TtT) and flTE| ) while (rrii) are by Eq. ( |3B| ) and all the fi(x) are normalized as 

( 1 dxf i {x) = l. (39) 
J o 

It is obvious that the r.h.s of Eq. (^), being integrated over x from to 1, just gives the 
average multiplicity of the final state protons for each reactions in the beam fragmentation 
region. 
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Three fragmentation functions can be determined if data on three different reactions are 



given. In addition to pp —> pX of Ref. [E3|, we use data on pCu — > pX and pAg — > pX 



provided by two experimental groups [T^J2^] in order to cover the wide x range. As Ref. 
f2?|| gives only the cross sections at fixed pr's, we assumed that they are proportional to the 
Pr-integrated cross section; 



da da 
x— = k E— 
dx dp 



(40) 

fixed pt 



Actually the constant k can be determined by requiring that the data of Ref. |27|] a t Pt = 0.3 
GeV/c coincide with the data of Ref. |1J in the range 0.3 < x < 0.6 where there exist data 
points from both groups. Such a procedure is possible because the x-dependence of a fixed 
Pr-cross section of [27j is really similar to xda/dx of |19[ as shown in Figs. ^| (c) and (d). 



The values of Kj £1X6 Kj — 1.30 and 1.20 [(GeV/c) 2 } for A = Cu and Ag, respectively. For 
information, the same procedure has been applied to pp — > pX with the result that k = 1.18 
[(GeV/c) 2 ]. See Fig. |6] (a). It is remarkable that the values of k are almost the same, i.e., 
independent of A in the three cases. 

Applying data on three different reactions to Eqs. (|38a| ) and ([38b]), one can determine 



the three fragmentation functions. However, the immediate application gives too large 
uncertainties for f'2{x) and fs(x) due to the propagation of the experimental errors. We 
accordingly assume that fs(x) ~ for x ^ 1/3 because a wounded nucleon having all the 
three interacted quarks loses most of incident momentum. The fragmentation functions thus 
obtained are shown as points with error bars in Fig. |5|, where they represent the values of 
the fragmentation functions including the flavor factors, Fi(x) = Xifi(x) at various x. 

Though the fragmentation functions have been determined at various x, we need Fi(x) 
for all x in order to calculate the nuclear stopping power in any reactions. In this case a 
meaningful x 2 fit f° r whole x region is impossible because of absence of experimental data 
at x < 0.2 and the large errors of data point at 0.2 < x < 0.3. Therefore, we assume 
appropriate functional forms and impose some physical conditions: Fi(x) is assumed to be 
a linear function of x, F2(x) and F%(x) are gaussian and exponential, respectively, with 
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conditions on the average fractional momenta (x)\ > (x) 2 > (x)^ and another condition 
on the flavor factors Ai > A2 > A3. Here, we do not use the pp — ► pX data point at 
x > 0.9 shown in Fig. |6|(a) because they are dominated by the diffractive contribution. 
The inequality for the average fractional momenta (x)i = J dxxfi(x) means that a wounded 
nucleon with the more interaction loses the more momentum. The probability of the flavor 
change in the incident nucleon is larger for the more interacted quarks giving the inequality 
for the flavor factors. The resultant fragmentation functions are 

F x (x) = 0.56 + 0.182, (41a) 
F 2 (x) = 1.4{exp(-3.54x 2 ) - exp(-3.54)}, (41b) 
F 3 (x) = 3.25exp(-6.52x). (41c) 

As shown in Fig. [5], Fi(x) is nearly constant while both F 2 (x) and Fs(x) decrease toward 
zero as x increases. The flavor factors are Ai = 0.65, A2 = 0.61 and A3 = 0.50 and the 
average fractional momenta are (x)i = 0.52, (x) 2 = 0.28 and (x)^ = 0.15. 

In Fig. ^|, the calculated proton spectra for pA — * pX reactions (A = p, Be, Cu, Ag, W 
and U) are compared with the experimental data [|TBJ [TP|] P7| . Our theoretical proton spec- 
tra estimated from Eq. (^) with above fragmentation functions Eq. ( f4"l~D and probabilities 
of quark absorption given in Table | reproduce not only the input data (p, Cu, Ag) but also 
the other data (Be, W, U). This results suggest the validity of CQM. 
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VI. RAPIDITY DISTRIBUTION OF PROTONS IN AA COLLISIONS 



The rapidity distributions of participant protons have recently been measured for S + S 
collisions at E lab = 200 AGeV by NA35 Collaboration j§ and also for central 208 Pb + 
208 Pb collision at Ei ab = 160 AGeV by NA44 Collaboration ||. In our model the fractional 
momentum distribution of protons can be calculated from Eq. (|38c|) with fi(x) given by 
Eq. (PD, (rrii) given in Table [TI| and the flavor factors A| given in Table [TTT]. Calculation 
of the factors is given in Appendix [FJ. The rapidity distribution can be obtained from 
the fractional momentum distribution ( |38cj ) using the relation 



where ttit is the transverse mass of the particle; m^ 2 = rn? + p\. 

The results for S + S collisions are shown in Fig. |7|. We have used the average values pt 
= 0.622, 0.595 and 0.45 GeV/c for central(veto and trigger) and peripheral collisions, 
respectively 0]. There is a good agreement between theoretical results and experimental 
data, particularly in the midrapidity region. The rise of the calculated spectra at very small 
yiab is due to the behavior of fi(x) at x ~ 1 which is rather ambiguous. There is a trend that 
the proton yield is slightly overestimated for the central triggers while it is underestimated 
for the peripheral collision. This may be due to a sharp cut of the impact parameter in our 
calculation. Our model can also reproduce well the preliminary data on central Pb + Pb 
collision || as shown in Fig. |^. Considerable yield of protons in the central rapidity region ( 
y lab ~ 3 ) is due to the dominance of (777.3 (m)/m) at the largest m. See Fig. §(b). 

In Fig. H the theoretical predictions is given for 197 Au + 197 Au collisions at ^/s = 200 
AGeV to be measured at BNL-RHIC. It is remarkable that the midrapidity region becomes 
nearly baryon free in contrast to Pb + Pb collisions at CERN-SPS(Fig. ^|). 

As a useful measure of nuclear stopping power, one usually uses the mean rapidity shift 
of the projectile proton from their original beam rapidity, (Ay) or the mean fractional 
momentum retained by the final state proton, (x) ||28|| . The larger the (Ay), the more 




(42) 



18 



stopping is present and hence the more baryons are produced in the central rapidity region. 
In Table [TV], (Ay) and (x) and the integrated yield (N p ) = J dy(dN/dy) calculated by our 
model are shown in comparison with experimental values. The rapidity shifts in S + 5* 
reactions for various triggers are in good agreement with the experimental values. The (Ay) 
in both central triggers are significantly larger than one in peripheral trigger, which implies 
the larger stopping in central collisions than in peripheral collisions. The rapidity shifts (Ay) 
in Au + Au collision at ^/s = 200 AGeV is much larger than that in Pb + Pb collision at 
17.4 AGeV, although (m) in two processes are similar to each other. This increase of (Ay) 
is due to the increase of (717,3). in spite of the increase of (Ay) at RHIC energy region, the 
baryon number density decreases at the central rapidity region because the beam rapidity 
also increases. The mean fractional momentum (x) in central collisions of heavy nuclei is 
smaller than that in central collisions of light nuclei by about 0.1. However, still more than 
20 % of the incident momentum is carried by the leading nucleons in central collisions of 
heavy nuclei. 

The mean proton multiplicity (N p ) calculated from our model for S + S collision is in 
good agreement with data within experimental errors. It is remarkable that (N p ) in central 
collisions of heavy nuclei is larger than the number of incident protons. This is due to the 
charge asymmetry of incident nuclei. Consider, for example, collisions of fictitious nuclei 
made of only neutrons. 
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VII. CONCLUSIONS AND DISCUSSIONS 



It has been shown that our theoretical results of proton spectra are in good agreement 
with experimental data for both the fractional momentum distribution in pA reactions and 
the rapidity distributions in central 32 S + 32 S collision at Ew^ = 200 AGeV and central 
208p b + 2os pb collision at = iQOAGeV (CERN-SPS). This result suggests strongly 

the validity of the constituent quark model. It is predicted that the central rapidity region 
in 197 Au+ 197 Au collision at RHIC energy region will be nearly baryon-free. In general the 
baryon number density in the central rapidity region increases with increasing mass number 
of colliding nuclei, whereas it decreases with increasing incident energy. During the course 
of this analysis, we have established a formula, Eq. ([35 ), for the probability of having 
three types of wounded nucleons in AB collisions and have determined the fragmentation 
functions, Eq. (f|l]), for those wounded nucleons. The probabilities given by Eqs. (|IT)|). 
(|18|) and (J35|), of having different quark interactions in NN, NA and AB collisions should 
be useful to evaluate various spectra, e.g., the transverse energy distribution, in nuclear 
collisions at quark level. 

It is worthwhile to compare our model with other models for nuclear stopping. For S + 
S reaction, the experimental data show a flat rapidity distribution of protons in agreement 
with our model and VENUS [0 while HIJING [TT|] shows stronger transparency. For Pb 
+ Pb collision, our model and RQMD JT(J| predict similar strong stopping || contrary to 
VENUS and HIJING which give weaker stopping. The different models lead to the different 
predictions for the nuclear stopping. It should be stressed that the proton spectra of AB 
collisions in our model result from fitting the pp and pA data through Eq. (|38|) . In general, 
any reasonable extrapolation from pp to AB via pA will give a similar result, a rather strong 
nuclear stopping. 

Our model can describe the proton spectra in pp, pA and AB collisions in a unified 
manner. The notable feature of the fundamental formula, Eq. ([38]), not shared with other 
models, is the factorization of x- and A-dependences. Equation fl3~8"D summarizes the crucial 
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feature of our dynamical assumptions that different types of wounded nucleons distinguished 
by the number of interacted quarks hadronize differently and independently. One has to 
notice that fragmentation functions which have been determined by the experimental data 
in this work should be eventually derived from theoretical considerations in future. Anyway, 
we would like to stress that the constituent quark model is so simple that it has only a 
few number of free adjustable parameters or functions; a™ el and fi(x) for i — 1, 2 and 3. 
Although most models have recently been constructed as complicated event generators, we 
feel it still worthwhile to pursue a simple phenomenological model which allows an analytical 
calculation and a simple interpretation of a result of data analysis. 
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APPENDIX A: THE PROBABILITIES OF HAVING j INTERACTED QUARKS 
IN n NUCLEON-NUCLEON COLLISIONS 



In NA collisions, we have considered that the incident nucleon collides with many nucle- 
ons of the target nucleus by multiple collision. By expanding the formula of total inelastic 



cross section Eq. (HQ, the probabilities of having n (nucleon-nucleon) collisions in an average 



collision with all possible impact parameters is written as 

Pna(u) = 4x / db( A \ [V A (b)} n [1 - V A (b)} A - n . (Al) 



To obtain the probabilities of projectile nucleon having j interacted quarks in n collisions, 
it is enough to expand only the second factor of Eq. (|A1|) [VA(b)] n , the probability of having 
exactly n collision, in terms of (6), the probability of finding i interacted quarks in one 
nucleon-nucleon collision, as follows; 

iv A (b)T = EE(Y) ( n 7 %vi 1 \b)nvi 2 \b)Y{v{ s \b)r- k - 1 . (A2) 



k=0 1=0 

r(i) 



The factor {V A (b)} on the r.h.s imply the probability of having k collision with the same 
V%\b) and can be expressed in terms of U^(k; 6), the probabilities of having j interacted 
quarks in k collisions while having i interacted quarks in each nucleon-nucleon collision. For 
k collisions with the same probability V^^b), it is given by 

{V£\b)} k = U^{k- b) + tf<g(fc; b) + 17<g(fc; 6), (A3) 



where 



U^(k ] b) = ?>^-V { 2\b))\ (A4) 
f/ ( ( 1 2 ) ) (fc;6) = 3(2«-2)(ivi 1) (6))", (A5) 
U$(k; b) = {(3" - 3) - 3(2" - 2)} (^(6))* ■ (A6) 

Here, Eq. ( |A4[ ) is the probability of one quark in the incident nucleon interacting repeatedly 
with any nucleon of target nucleus and then having one interacted quark in k collisions. Eq. 
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implies the probability of one quark interacting with the target nucleon, also another 
quark interacting in other collision, then having two interacted quarks. The probability 
of three quarks interacting in each collision is appeared in Eq. ( |A6| ) . When the projectile 
nucleon interacts with the target nucleon simultaneously in one collision, there are two 
possibilities of having two interacted quarks and three interacted quarks in / collision as 
follows; 

{V^b)} 1 = U$(hb) + U$frb), (A7) 

where 

r(2)/,. l\ o ( l \A 2 )( 



U^(l;b)=3^-Vr(b)j , (A8) 



U$(l;b) = {3»-3) (-Vf (6)) . (A9) 



In case where all the three incident quarks are participating in one N — N collision, the 
projectile nucleon has necessarily three interacted quarks in m N — N collisions; 

{vP(b)r = U^(m;b), (A10) 



Substituting Eqs.flA3|), flA7|) and ([AToD for Eq. Q, we obtain 



[VA(b)] n = EE Q ( n l k ) {*/$(*; b) + u$(h, b) + u$(k-, &)} 

x {U$(l- b) + U$(l; b)}{U$(n — k — I; 6)}, (All) 

In Eq. ( |A11| ), all terms including j = 3 contribute the probability having three interacted 
quarks in n collisions and the cross terms are divided into the probabilities of having two 
or three interacted quarks with suitable weight. Thus one can write [VA(6)] ra in terms of 
Ujp A (n;b), the probability of having j interacted quarks in n nucleon-nucleon collision in 
the following form, 

[V A (b)] n = uS\(n; b) + uP A (n; b) + U% ] A {n- b) (A12) 

where 
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uS\fab) = U$(n-,b), 



U^(n ] b) = U%(n ] b) + Ug>(n-,b) 



(i) 

n-l 

+ E 

k=l 



n 



U^Uk; b)U$(n -k;b) + ^Uft(k; b)U%{n - k; b) 



3^(i) 



r(2) 



(2), 



r(2), 



'(2) 



'(2) 



(A13) 

(A14) 
(A15) 



From Eqs.( |A12] ) and ( |A1| ), therefore, we obtain the formulas for probability function of 
quark absorption in average n collision . 



APPENDIX B: CALCULATION OF THE FLAVOR FACTORS IN AA 

COLLISIONS 



The flavor factors in AA collisions is the weighted average of the probabilities that a 
proton is produced from an incident proton or from an incident neutron. If one neglects the 
effect of baryon-antibaryon pair production and assumes that an incident nucleon fragments 
into either p, n or A one has 



A 



(A) 



(Bl) 



where rj is the probability that an incident nucleon is converted into A after interaction 
and Z/A is the proportion of the proton over the atomic mass. For a charge symmetric 
system(ZA4 = 1/2), it reduces to 

1 



\[ A) = \ { 2 A) 



A 



(A) 



(B2) 



The numerical value of rj is estimated to be 0.17 by using the experimental data on central 
S + S collisions @,|29 . 



24 



REFERENCES 

[1] L. D. Landau, Izv. Akad. Nauk. SSSR, Ser. Fiz. 17, 51(1953). 
[2] J. D. Bjorken, Phys. Rev. D 27, 140(1983). 

[3] J. Barrette et al., Phys. Rev. Lett. 64, 1219(1992); ibidem, 70, 2996(1993). 
[4] J. Bachler et al., Phys. Rev. Lett. 72, 1419(1994). 

[5] S. Esumi, talk presented at International Symposium on Quantum Interferometry Stud- 
ies in High Energy Nuclear Collision, Hiroshima, Japan, 1995; Genshikaku kenkyu 40, 
159(1996). 

[6] R. J. Glauber, in Lectures in Theoretical Physics, edited by W. E. Brittin and L. G. 
Dunham(Interscience, New York, 1959), Vol. I, p. 315. 

[7] K. Werner, Phys. Rep. 232, 87(1993). 

[8] See for a review, C. Y. Wong, Introduction to high-energy heavy-ion collisions, World 
Scientific, Singapore, 1994. 

[9] M. Gyulassy, Nucl. Phys. A590, 431c(1995). 

[10] H. Sorge, H.Stocker and W. Greiner, Nucl. Phys. A498, 567c(1989); ibidem, A566, 
633c(1994). 

[11] X. N. Wang and M. Gyulassy, Phys. Rev. D 44, 3501(1991); ibidem, 45, 844(1992). 
[12] V. V. Anisovich, Yu. M. Shabelski, and V. M. Shekhter, Nucl. Phys. B133, 477(1978). 
[13] N. N. Nikolaev and S. Porkorski, Phys. Lett. B80, 290(1979). 
[14] A. Dar and F. Takagi, Phys. Rev. Lett. 44, 768(1980). 
[15] A. Bialas and E. Bialas, Phys. Rev. D 20, 2854(1979). 
[16] F. Takagi, Prog. Theor. Phys. 65, 1350(1981). 

25 



[17] F. Takagi, Phys. Rev. D 27, 1461(1983). 

[18] A. E. Brenner et al, Phys. Rev. D 26, 1497(1982). 

[19] R. Bailey et al, Z. Phys. C 29, 1(1985). 

[20] G. G. Simon et al, Nucl. Phys. A333, 381(1980). 

[21] Review of Particle Properties, K. Hikasa et al, Particle Data Group, Phys. Rev. D 45, 
SI (1992), 

[22] A. Bialas, W. Czyz and L. Lesniak, Phys. Rev. D 25, 2328(1982). 

[23] A. Bohr and B. R. Mottelson, Nuclear structure, Vol. 1: Single- Particle Motion, W. A. 
Benjamin, Inc, New York, 1969. 

[24] S. P. Denisov et al, Nucl. Phys. B61, 62(1973). 

[25] A. S.Carroll et al, Phys. Lett. B80, 319(1979). 

[26] E. O. Abdrahmanov et al, Z. Phys. C 5, 1(1980). 

[27] D. S. Barton et al, Phys. Rev. D 27, 2580(1983). 

[28] See for example S. Date, M.Gyulassy, and H. Sumiyoshi, Phys. Rev. D 32, 619(1985). 
[29] J. Bartke et al, Z. Phys. C 48, 191(1990). 



26 



FIGURES 

FIG. 1. Three possible interactions of the incident nucleon colliding with the target, nucleon 
or nucleus, in CQM; one(a), two(b) and all the three(c) incident quarks are interacted with the 
target. 

FIG. 2. The A dependence of the probabilities that one, two or three quarks in an incident 
nucleon interact in pA collision. 

FIG. 3. The probability distribution Pab(^) of having m wounded nucleons for different trig- 
gers in (a) 32 S + 32 S" and (b) 208 Pb + 208 Pb. Here we have used of^ = 4.32 mb. 

FIG. 4. The proportion of three types of wounded nucleons rtij to total wounded nucleons m 
in no-bias event of (a) 32 S + 32 S and (b) 208 Pb + 20S Pb collisions. 

FIG. 5. The fragmentation functions Fi(x) obtained from the fractional momentum distribution 
of protons in pp, pCu and pAg reactions. 

FIG. 6. The proton spectra in pA collisions for A = p(a), Be(b), Cu(c), Ag(d), W(e) and U(f); 
pp, pCu and pAg data are used as inputs. Curves show the model result while data points are 



taken from [18|( py-integrated pp data at 100 GeV/c; triangles), JT9|( p^-integrated pA data at 



100 GeV/c; solid circles) and |27| ( p-p-fixed pp and pA data at 120 GeV/c; open squares). 

FIG. 7. The rapidity distribution of participant protons in 32 S + 32 S collision at Ei^ = 200 
AGeV The solid, broken and dotted lines correspond, respectively, to Veto, Et and periph- 
eral triggers. The corresponding experimental data are shown the circles, triangles and squares, 
respectively. 

FIG. 8. The proton rapidity distribution in central (15% trigger) 208 Pb + 208 Pb collision at 
Eiab = 160 AGeV . Solid line stands for theoretical result which includes the contribution from A 
decay. Experimental data are taken from Ref. ||. 
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FIG. 9. Theoretical prediction for proton rapidity distribution in central (15% trigger) 
197 Au + 197 Au collision at y/s = 200 AGeV. Contribution from A decay is included. It amounts to 
some 20% of the total yield. 
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TABLES 

TABLE I. Comparison of a calculated total inelastic cross section of^, with the experimental 
values from Ref. [|19| and |25| |, the probabilities having one, two or all the three interacted 
quarks in pA collisions and their average value < i > . 



A 


°!nel\cal(mb) 


<? P inel\exp{mb) 


P (i) 
r pA 


p(2) 
r pA 


p(3) 
r pA 


< i > 


P 


30 


31.3±1.2 


0.81 


0.17 


0.02 


1.21 


9 Be 


188 


176±2 


0.61 


0.30 


0.09 


1.48 


32 S 


493 




0.47 


0.34 


0.19 


1.72 




785 


767±8 


0.39 


0.34 


0.27 


1.88 


108 Ag 


1118 


1097±12 


0.33 


0.33 


0.34 


2.01 


189 W 


1584 


1540±16 


0.27 


0.31 


0.42 


2.15 


208p 6 


1724 


1752±53 


0.26 


0.30 


0.44 


2.18 


238 jj 


1880 


1860±20 


0.25 


0.29 


0.46 


2.21 
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TABLE II. Comparison of a calculated total inelastic cross section cr^A with the experimental 
values from Ref. Q and empirical formula, the average number of wounded nucleons (m) and the 
average values for three type of wounded nucleon {mi), (m?), (ins) for different trigger in AA 
collision. Here "no-bias" implies that the impact parameter integration is carried out from zero to 
infinity. 



AA(s/s(AGeV)) 


trigger 




<r££l\ex P (mb) 


(m) 


(mi) 


(m 2 ) 


(m 3 ) 


S+S(19.4) 


no-bias 


2173 


1700(1740*) 


9.4 


5.6 


2.9 


0.9 




Veto(2%) 


43 


34 


27.2 


9.8 


11.1 


6.3 




E T (ll%) 


239 


190 


25.3 


10.6 


9.9 


4.8 




periphe. 


1422 


1000 


4.0 


3.1 


0.8 


0.1 


Pb+Pb(17.4) 


no-bias 


7660 


7630* 


77.0 


32.7 


26.6 


17.8 




cent. (15%) 


1149 




196.0 


39.5 


75.5 


81.0 


Au+Au(200) 


no-bias 


7414(7179 t ) 


7325* 


80.2 


25.8 


25.8 


28.6 




cent. (15%) 


1112 




192.7 


18.2 


54.0 


120.5 



Values of empirical formula in Ref. [26|, see text. ' Value for (T^ei = ^-32 



TABLE III. Flavor factors \[ A) calculated by Eqs. PD and in Appendix 



A 



(A) 
1 



A 



(A) 



A 



(A) 



s 
Pb 

Au 



0.415 
0.365 
0.368 



0.415 
0.374 
0.376 



0.415 
0.397 
0.398 
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TABLE IV. The mean multiplicities (N p ) and the mean rapidity shift (Ay) for y < y cm of the 
final state proton in our theoretical results are compared with the experimental values in Ref. j| 
for S + S reaction at Ei a b = 200 AGeV. (x) is the mean fractional momentum retained by the final 
state proton. 



AA(y/s(AGeV)) 


trigger 


{Np/ 1 exp 




(Ay p )\ exp 




< x > \ cat 


S+S(19.4) 


Veto(2%) 


12.8±1.4 


14.4 


1.58±0.15 


1.52 


0.34 




E T {11%) 


10.3±1.4 


13.3 


1.58±0.15 


1.46 


0.36 




periphe. 


3.1±0.8 


2.0 


1.00±0.15 


1.16 


0.47 


Pb+Pb(17.4) 


cent.(15%) 




96.4 




1.67 


0.27 


Au+Au(200) 


cent.(15%) 




96.6 




2.22 


0.22 
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(a) 32 S + 32 S 
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(a) 32 S + 32 S 



1 .0 I i i i i I i i i i I i i i i I i i i i I i i i i I i i i 



0.8 



0.6 



0.4 



0.2 



0.0 



<m i (m)>/m 

— — - <m 2 (m)>/m 
<m 3 (m)>/m 







10 15 20 25 30 

m 

Fig.4(a) 



1.0 



0.8 



0.6 



0.4 



0.2 



0.0 



(b) 208 Pb + 208 Pb 
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32 S + 32 S protons (E = 200 AGeV) 
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